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We demonstrate that the requirement of galilean invariance determines the choice of H function 
for a wide class of entropic lattice Boltzmann models for the incompressible Navier-Stokes equations. 
The required H function has the form of the Burg entropy for D — 2, and of a Tsallis entropy with 
q = 1 — for D > 2, where D is the number of spatial dimensions. We use this observation to 
construct a fully explicit, unconditionally stable, galilean invariant, lattice-Boltzmann model for the 
incompressible Navier-Stokes equations, for which attainable Reynolds number is limited only by 
grid resolution. 



INTRODUCTION 

Lattice Boltzmann models of fluids (T], [| evolve a 
single-particle distribution function in discrete time steps 
on a regular spatial lattice, with a discrete velocity space 
comprised of the lattice vectors themselves. The single- 
particle distribution corresponding to lattice vector Cj at 
lattice position x and time step t is denoted by iVj(x, t). 
The simplest variety of lattice Boltzmann models employ 
a BGK operator ||, so that their evolution equation is 

Ni(x + c u t+ At) = Nifa t) + - [iV^ q (x, t) - 7V 4 (x, t)] 



for i = 1, . . . , b. Here b is the coordination number of the 
lattice, iVj Cq (x, t) is a specified equilibrium distribution 
function that depends only on the values of the conserved 
quantities at a site, and r is a characteristic collisional 
relaxation time. Using the Chapman-Enskog analysis it 
is possible to show that the mass and momentum mo- 
ments of the distribution function will obey the Navier- 
Stokes equations for certain choices of equilibrium distri- 
bution J|. 

The viscosity that appears in the Navier-Stokes equa- 
tions obtained from these models is proportional to r — \ . 
To lower viscosity and thereby increase Reynolds num- 
ber, practitioners often over-relax the collision operator 
by using values of r in the range (£, 1]. For sufficiently 
small r, however, the method loses numerical stability, 
and this consideration limits the lowest Reynolds num- 
bers attainable. 



In an effort to understand these instabilities, there has 
been much recent interest in entropic lattice Boltzmann 
models ||, ||, |(| . These models are motivated by the fact 
that the loss of stability is due to the absence of an H the- 
orem. Numerical instabilities evolve in ways that would 
be precluded by the existence of a Lyapunov function. 
The idea behind entropic lattice Boltzmann models is to 
specify an H function, rather than just the form of the 
equilibrium. Of course, the equilibrium distribution will 
be that which extremizes the H function. The evolution 
will be required never to decrease H, yielding a rigor- 
ous discrete-time ff-theorem; this is to be distinguished 
from other discrete models of fluid dynamics for which 
an //-theorem may be demonstrated only in the limit of 
vanishing time step [Q. 

To ensure that collisions never decrease H, the collision 
time r is made a function of the incoming state by solv- 
ing for the smallest value r min < 1 that does not increase 
H. The value then used is r = r min /«; where < K < 1. 
It has been shown that the expression for the viscosity 
obtained by the Chapman-Enskog analysis will approach 
zero as k approaches unity |i[ ^, g|. Thus, the entropic 
lattice Boltzmann methodology allows for arbitrarily low 
viscosity together with a rigorous discrete-time H theo- 
rem, and thus absolute stability. The upper limit to the 
Reynolds numbers attainable by the model is therefore 
determined by loss of resolution of the smallest eddies, 
rather than by loss of stability |^|, ^| . 

In a recent review of the subject, Succi, Karlin and 
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Chen [[10| have pointed out that entropic lattice Boltz- 
mann models have three important desiderata: Galilean 
invariance, non-negativity of the distribution function, 
and ease of determining the local equilibrium distribu- 
tion at each site at each timestep. 

In this paper, we shall construct entropic lattice Boltz- 
mann models for the incompressible Navier-Stokes equa- 
tions which are gallilean invariant to second order in 
the Mach number expansion of the distribution func- 
tion (quasi-perfect in the terminology of ]l0||). We shall 
show that the requirement of galilean invariance makes 
the choice of H function unique. We shall show that the 
required function has the form of the Burg entropy |Tl] ] 
in two dimensions, and the Tsallis entropy in higher di- 
mensions. While the analogous problem for the com- 
pressible Navier-Stokes equations is difficult and remains 
outstanding, the purpose of this paper is to point out 
that the incompressible case is nontrivial and interesting 
in its own right. 

EQUILIBRIUM DISTRIBUTION 

We consider a Bravais lattice of coordination number 
b in D dimensions. We denote the lattice vectors by Ci 
where i = 1, . . . , b, and their magnitudes by c = |cj|. We 
demand that the lattice symmetry group be sufficiently 
large that the only fourth-rank tensors that are invari- 
ant under its group action are isotropic. The mass and 
momentum densities are given by 

b 

1=1 

and 

b 

pu = y^mCjiVi, (2) 

i=l 

where m is the particle mass, and u is the hydrodynamic 
velocity D-vector. These D + 1 quantities must be con- 
served in collisions. 

If we regard the JVj, for i = 1, . . . , b, as coordinates in 
a ^-dimensional space, the conservation laws ([!]) and (||) 
restrict the collision outcomes to a b — (D + 1) dimen- 
sional subspace. Since the conserved quantities are linear 
functions of the iVj's, the nonnegativity requirement 

Ni > (3) 

is satisfied within a compact polytope whose faces are 
given by the b equations Ni = for i = 1, . . . , b. We 
assume that the ii-function is of trace form 

b 

i=i 



where h'(x) > for x > 0. If lim^^o h'(x) = oo, 
then the normal derivative of H goes to negative infin- 
ity on the polytope boundary, enforcing the nonnegativ- 
ity constraint, Eq. (||). The purpose of this paper is to 
demonstrate that the requirement of galilean invariance 
uniquely determines the choice of function h(x). 

The equilibrium distribution function may be found 
by extremizing H with respect to the Ni, subject to the 
constraints, Eqs. ([!]) and (0), 

ol\i \ m m ) 

where fi/m and (3/m are Lagrange multipliers. We 
quickly find 

o = ti (Nd-n-p-a, 

and so 

A^ = 0( M + /3- Ci ), (4) 

where the function <f> is the inverse function of h! . The 
constants fi and (3 are determined by Eqs. (Q) and (Q), 
though it is generally impossible to find an exact analytic 
expression for them in terms of the conserved quantities 
p and pu; rather one must solve for them numerically or 
perform a Taylor expansion in Mach number. We adopt 
the latter approach below. 

GALILEAN INVARIANCE 

We seek to Taylor expansion the equilibrium distri- 
bution in Mach number because (i) we can do so ana- 
lytically, (ii) only the first two terms of that expansion 
determine the form of the incompressible Navier-Stokes 
equations, and (iii) that expansion is a useful initial guess 
for any numerical solution. From general symmetry ar- 
guments it is clear that (3 will be proportional to the hy- 
drodynamic velocity u, so that we may begin our Mach 
number expansion by expanding Eq. (Q) for small (3. We 
get 

N- q = 0(m) + ■ ^ + ^"(/W : c iCi + ■■■ 

Inserting this into Eqs. (0) and (||), and using general 
properties of the Bravais lattice, we find 

p = m b^{n) + —4>"{ix)P 2 + --- 

and 

pu= _^_0'( At )/3 + ..., 

where the ellipses denote third or higher order in Mach 
number. Inverting this perturbatively we find that, to 
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second order in Mach number, the Lagrange multipliers Inserting these into Eq. (||), we obtain the equilibrium 
are given by distribution, 

„_ x 

2 « 2 WW? 
where x = h! (^5), and by 

n - E - 

f3 = 2ifc_ u ^ . 

c 2 4>' (x) 
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mb 
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1 + -rrCi ■ U 



D 2 <t> (x)(f>" (x) 



2c 4 [^{x)] 2 V s - D 
Now it is well known that a Chapman-Enskog analysis based on the equilibrium distribution 
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(G) 



will give rise to the incompressible Navier-Stokes equa- 
tions 



and 



du 

~dt 



V-u = 



gu ■ Vu = — VP 

P 



vV 2 n. 



Comparing Eqs. (||) and (|J), we identify 

D \ <f> (x) 4>" (x) 
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The factor g destroys the form of the convective deriva- 
tive, and hence breaks galilean invariance. To recover 
galilean invariance we demand that g = 1 , and this yields 
the second-order nonlinear differential equation 



W(x)X 



The general solution to this equation is of the form 

(j>{x) =C D/2 {x-aC) p , 

where C and a are arbitrary constants, and P is to be 
determined. We quickly find that P must be either or 
—D/2. Since a constant (j> would not yield a well defined 
hi , we see that we must have <f>(x) — C D ^ 2 (x ~ aC)~ D l 2 , 
whence h'(x) = C(a + x~ 2 / D ), and this integrates to give 



h(x) 
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if D = 2 
if D ^ 2, 



(7) 



where h is constant. In fact, the only effect of nonzero 
ho is to introduce an additive constant to H, and the 
only effect of nonunity C is to scale H by a constant 
factor. In other words, h(x) is uniquely specified only 
to within additive and multiplicative constants. With 
this understanding, we may say that the requirement of 
galilean invariance has uniquely specified the choice of 
H. 

In passing, we note that linr^o h r (x) = 00. Thus the 
nonnegativity constraint will be enforced by the dynam- 
ics. 

Finally, we write the global Lyapunov function TL = 
J2*H by summing h(Ni(x,t)) over the lattice. Since 
the total mass is conserved we have complete freedom 
to choose a, and so to within additive and multiplicative 
constants TL may be written 

„ m f£x£>[^(^)] fori) = 2 

" \ £ x £ 4 for D jt 2, 

for appropriate choices of a and C . This has the form of 
a Burg entropy JTTJ for D = 2, and a Tsallis entropy 
with parameter 

2 

for D ^ 2. We note that D < 2 corresponds to q < 0, 
and D > 2 corresponds to q > 0. It is interesting that it 
is only in the infinite-dimensional limit, D — > 00, where 
the set of velocities becomes infinite, that q — > I and we 
recover the Boltzmann-Gibbs entropy |l3| . 

The appearance of the Burg and Tsallis entropies in 
this context is fascinating. In a footnote of their re- 
cent review, Succi, Karlin and Chen noted that the 
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entropy that gave rise to the above-mentioned solvable 
model for a compressible fluid was related to the Tsal- 
lis entropy with q = 3/2, so there may be more than one 
connection with Tsallis thermostatistics |Q lurking here. 
There are precious few situations in which the origins of 
Tsallis thermostatistics can be traced analytically to an 
underlying microscopic model, as we have done here. 
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CONCLUSIONS 

We have presented galilean- invariant, entropic lattice 
Boltzmann models for the incompressible Navier-Stokes 
equations. We expect that these models will be useful for 
the simulation of two- and three-dimensional turbulence. 
As noted by Succi, Karlin and Chen pToj ] , the problem of 
finding perfect models for lattice models of the compress- 
ible Navier-Stokes equations is much more difficult and 
may well be impossible. We found it interesting that the 
simpler problem, for incompressible fluids, is itself very 
nontrivial and interesting. In particular, the appearance 
of the Burg and Tsallis entropies for the H function is 
surprising. These entropies have heretofore been asso- 
ciated with long-range interactions, long-time memory 
or a fractal space-time structure. This work indicates 
that they may also be relevant to models with discretized 
space-time, and this surely warrants future study. 



Acknowlegements 

BMB was supported in part by the U.S. Air Force Of- 
fice of Scientific Research under grant number F49620- 
01-1-0385, and in part by MesoSoft Corporation. He 
performed a portion of this work while at the Centre 



[1 
[2 
[3; 
[4 
[5] 
[6] 

ir 

[8 

[9 
[10 

[11 
[12 

[13 



S. Succi, The Lattice Boltzmann Equation - For Fluid 
Dynamics and Beyond (Oxford University Press, 2001). 
R. Benzi, S. Succi, and M. Vergassola, Phys. Reports 222 
(1992). 

Y. H. Qian, D. d'Humieres, and P. Lallemand, Europhys. 
Lett. 17, 479 (1992). 

I.V. Karlin, A.N. Gorban, S. Succi, and V. Boffi, Phys. 
Rev. Lett. 81, 6 (1998). 

I.V. Karlin, A. Ferrante, and H.C. Ottinger, Europhys. 
Lett. 47, 187 (1999). 

B. M. Boghosian, J. Yepez, P. V. Coveney, and A. Wag- 
ner, Proc. R. Soc. (London) A. 457, 717 (2001). 

C. Marsh and P. V. Coveney, J. Phys. A Math. Gen. 31, 
6561 (1998). 

S. Ansumali, I.V. Karlin, Phys. Rev. E 62, 7999 (2000). 
S. Ansumali, I.V. Karlin, Phys. Rev. E 65, 056312 (2002). 
S. Succi, I.V. Karlin, H. Chen, Rev. Mod. Phys. (2002), 
to appear. 

A.N. Gorban, I.V. Karlin, Phys. Rev. E. accep ted for 



publication (2002), Preprint, |cond-mat/0205511 
C. Tsallis, Nonextensive Statistical Mechanics and Its 
Application (Springer Verlag, 2001), edited by S. Abe 
and Y. Okamoto. 

This limit must be taken with great caution, since the 
sound speed, c s — c/V^D, also vanishes in this limit, 
threatening the validity of the expansion in Mach num- 
ber M — u/c s ; we could, for example, take u/c ~ D~ 3 ^ 4 
so that M — > even as D — > 00. 



